Effect of Microwaves on the Current-Phase-Relation in SNS Josephson Junctions 
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We investigate the current-phase-relation (CPR) of long diffusive superconductor-normal metal- 
superconductor (SNS) Josephson junctions in thermodynamic equilibrium and under microwave 
irradiation. While in equilibrium good agreement with the predictions of quasiclassical theory is 
found, we observe that the shape of the CPR can be strongly affected by microwave irradiation. 
Close to a Josephson-phase difference p ~ tt, the supercurrent can be strongly suppressed when 
increasing the rf-power. Our results can be understood in terms of microwave excitation of low- 
lying Andreev bound states across the mini-gap in the junction. In the frequency interval studied, 
this mechanism becomes important, when the mini-gap closes at p ~ n. 

PACS numbers: 74.45.+C, 74.50.+r, 85.25.Cp 



The fundamental mechanism of Cooper-pair trans- 
fer across superconductor-normal metal-superconductor 
(SNS) Josephson junctions relies on the existence of An- 
dreev bound states (ABS) in the normal metal [H. The 
energies e of the ABS depend on the phase difference p 
between the two superconductors. Diffusive normal met- 
als exhibit a continuous spectrum of ABS, described by 
a complex supercurrent spectral density js(e,ip) [3, 0j|. 
Assuming perfectly transparent interfaces at a distance 
d, js(e, f) is determined by two parameters, namely the 
energy gap A of the superconductor, and the Thouless 
energy e t h = hD/d 2 , where D is the diffusion constant. 
The ABS spectrum exhibits a minigap e s (ip), which closes 
as tp approaches tt. The total supercurrent is given by Q 

H<P) = ^rl lm JS (e,p)[l-2f(e)}de, (1) 

where is the resistance of the normal conductor and in 
thermal equilibrium /(e) is the Fermi distribution func- 
tion determining the occupation probability of the ABS 
at temperature T. As opposed to the case of conven- 
tional tunneling junctions, for ideal NS interfaces higher 
order scattering processes become important. These 
processes correspond to the phase-coherent transfer of 
multiple Cooper pairs across the junction, and mani- 
fest themselves in a non-trivial Fourier representation 

I s (T,ip) = Y, j I c\ T ) sin U l P) of E q- L In thermody- 
namic equilibrium, the higher harmonics Iq^ are rapidly 
suppressed for temperatures larger than eth/^B- 

Attempts to detect higher harmonics in the CPR of 
diffusive SNS-junctions by measuring the magnetic re- 
sponse of loops with an embedded junction date back to 
the 70 's [H, but failed because of the contradictory re- 
quirements of large eth (very low temperatures and short 
junctions) and low Iq (higher T and long junctions) to 
avoid hysteretic switching of the loop Q . As an alterna- 
tive method, the detection of subharmonic Shapiro steps 
was suggested [1, 0]- These experiments were success- 
ful, but resulted in a very surprising non-monotonic In- 
dependence of the higher order coefficients Iq . As op- 
posed to the monotonic exponential suppression with T 



expected from theory [3( , the subharmonic Shapiro steps 
were observed most clearly at T > T t h- At the highest 
temperatures even I@ < ify and 1^ < I c 3 ^ were found 
0] ■ These findings have so far remained in contradic- 
tion with the equilibrium quasi-classical theory and were 
explained by the generation of non-equilibrium by the 
dc-bias required to records the IV characteristics. 

In this Letter, we present an investigation of the cur- 
rent phase relation in long diffusive SNS junctions using 
Waldram's method [!, 0| with a micro-Hall sensor Q to 
detect the magnetic flux. In thermodynamic equilibrium, 
we find quantitative agreement with the quasi-classical 
theory. In particular, the equilibrium CPR was found to 
be entirely sinusoidal at higher temperatures. When the 
samples are irradiated with microwaves with a photon 
energy comparable to e g (y), we could induce a signifi- 
cant contribution of higher harmonics in the CPR with 
an anomalous T-dcpcndcncc. Our observations can be 
understood in terms of a non-equilibrium occupation of 
Andrccv-levels in the microwave field, and are important 
for the control of Andreev-quantum bits Q . 

To investigate the supercurrent response of our SNS- 
junctions to an applied phase difference, we integrated 
the junctions into superconducting loops. The gauge in- 
variant phase difference p across the junction is then re- 
lated to the total magnetic flux <£> via p> = — 27r<J?/$o, 
where $0 = h/2e is the magnetic flux quantum [l 01 ] . 
The samples were patterned on top of the active area 
of micron-sized Hall-crosses (see insets in Fig. [IJ which 
have been structured into a GaAs/AlGaAs semiconduc- 
tor heterostructure containing a 2-dimensional electron 
gas 190 nm below the surface. With a mean free path of 
I = 9 /im at 4 K and an active area of 10 x 10 /jm 2 these 
Hall probes work close to the ballistic regime. Depend- 
ing on the bath temperature, which sets an upper limit 
to the applied current through the Hall probe, I p , wc 
obtain a supercurrent sensitivity between 0.1 and 2/xA, 
corresponding to a flux sensitivity of 10~ 3 to 2 • 10 -2 $o- 
The fabrication process makes use of e-beam lithography 
and a shadow evaporation technique based on a ther- 
mostable trilayer mask system. We use a sacrificial layer 
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FIG. 1: (color online) The measured Hall voltage Uh as a 
function of the loop current 7 cxt . For low temperatures (/3 > 
1) the effective flux shows hysteretic behavior, hence only the 
rising branch of the CPR is accessible. Upper inset: SEM 
image of the Nb loop with an edge length of 10 /im and an 
inductance L ~ 40 pH fabricated on the active area of a micro- 
Hall-cross. Lower inset: close-up of the actual SNS junction. 
The thickness of the Ag (Nb) layer is 40 nm (150 nm). While 
the spacing between the Nb electrodes is d — 495 nm, the 
effective length of the Ag bridge is d e fi = 640 nm with a width 
of 210 nm. 
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FIG. 2: (color online) CPR for low temperatures down to 
ksT ~ e t h- For T < 500 mK a clear deviation from the si- 
nusoidal Josephson relation is observed. The solid lines cor- 
respond to the predictions of the quasiclassical theory, with 
all parameters determined independently. Inset: The critical 
current as a function of temperature. From the fit of the data 
to the quasiclassical theory (solid line) we determine the value 
D = 0.014 m 2 /s of the diffusion constant in the Ag bridge. 



of polyether sulfone (PES) as described in the work of 
Dubos et al. [ll[ and a mask layer of silicon nitride (SiN) 
[13]. The insets to Fig. [1] show scanning electron micro- 
graphs of the loop on top of the Hall cross, and a close-up 
of the junction, respectively. 

The magnetic flux is controlled by driving an external 
current J ex t through the loop via current leads close to 
the SNS contact (Fig. [T] upper inset). This has the ad- 
vantage that both the external flux LI ey ± and the flux 
Lis induced by the junction are detected with the same 
filling factor a, which greatly facilitates the calibration 
of the measured Hall voltage J7h in terms of $o : 



the Hall cross, a' = aAnej! ] 



pj 



1.602- 



-19 



2,26 ■ 



10 15 m -2 , and L the geometric inductance of the loop. 
Depending on the match between the loop area and A, 
a-values between 0.3 and 0.4 can be achieved, limited by 
the finite distance between the 2DEG and the loop. 

Figure [T] shows the Hall voltage Uu(Iext) of a sample 
with d = 495 nm for different temperatures. The mag- 
netic response of the SNS loop is superimposed on a linear 
background caused by I C xt , the slope of which can be used 
to determine L. Below a certain temperature, when the 
parameter /?l = 2irLIc(T)/§ exceeds 1, the response is 
hysteretic with distinct jumps when the critical current 
of the junction is exceeded. In the limit ft,>l, the step 
height corresponds to one flux quantum and can thus be 
used to determine a'. A plot of Is = iext — a'Un/L vs. 



tp = 27to/(7h/ < I > o represents the desired CPR. 

The Thouless energy e t h of our junctions is determined 
by the effective length d c ff > d of the normal conducting 
barrier, which accounts for the fact that the Andrcev 
reflections at the SN interface occur randomly within the 
overlap length s between Nb and Ag. As in earlier work 
[13 ]. we have chosen d c s = d + s. The gap parameter 
A = 1.3 meV of our Nb-films has been extracted from its 
measured T c ~ 8.5 K. Fig. [2] shows the measured CPR 
for several temperatures below 1 K. While at 1 K the 
CPR is still a perfect sine, at lower temperatures clear 
deviations from the sinusoidal shape are observed. The 
inset to Fig. [2 shows Iq(T), extracted from the maximum 



$ = a$ H = a'U n = L (J ext - I s (27r*/$ )) ■ (2) supercurrent 



Here, $h is the flux enclosed by the active area A of 
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which agrees very well with the quasi- 
classical theory [3j, [l5j] using D = 0.014 m 2 /s as the only 
fitting parameter (solid line in inset). This value of D 
corresponds to a eth = 24 /ieV for this sample (d e e = 
640 nm). With all parameters fixed, we are now prepared 
for a quantitative comparison of the measured CPR with 
the theoretical predictions 0], represented by the solid 
lines in Fig. [5J Within our experimental resolution the 
agreement is excellent [ToT ]. In particular, the result for 
our lowest accessible temperature of 300 mK = 1.1 eth/ks 
almost coincides with the theoretical predictions for this 
value of A/eth = 58 in the limit T — > 0, where the degree 
of anharmonicity is largest. 

We now turn to the main subject of the present paper 
- the effect of microwave irradiation on the CPR in the 
frequency range between 3 and 12 GHz. The rf-power 
is applied at room temperature and reduced by several 
attenuators at 4K, IK and 100 mK (in total — 60 dB in- 
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FIG. 3: (color online) a) Measured CPR under irradiation 
with microwaves at T = 2.8 K and / = 7.0 GHz. b) Measured 
CPR for different temperatures at / = 11.8 GHz. Solid lines 
represent the best fit to our model. Inset: Similar data at 
7.4 GHz (near a cavity resonance). The solid line is a Fourier- 
expansion with the three coefficients listed. 



eluding cable losses) to ensure proper thermalization of 
the rf-cable, which is terminated with a small loop in the 
vicinity of the sample. Figure [3^ shows the CPR mea- 
sured at 2.85 K for different levels of the externally ap- 
plied power P. These data are taken away from a cavity 
resonance of the sample chamber. At this temperature, 
the equilibrium CPR for P = is nearly sinusoidal and 
the critical current is small enough < 1) to access 
the full CPR. It is clearly seen that with increasing power 
the maximum supercurrent is reached for lower tp- values, 
tpmax < tt/2, in contrast to the behavior at low tempera- 
tures (see Fig. [2]). The effect is most prominent for phase 
differences close to n, where a pronounced suppression of 
the supercurrent occurs. For tp < tp max , the supercurrent 
even appears to be slightly enhanced compared to P = 0. 
At certain frequencies, which we attribute to resonances 
of our sample chamber, we found that Iq and 1^ can 
even be comparable (see inset in Fig.[3Ji). This can even 
lead to an inversion of the supercurrent direction close 
to tp < ir. In Fig.[3jD the temperature dependence of the 
effect is investigated at P = 28 dBm. The deviation from 
the sinusoidal form is moderate at 2.3 K and becomes 
more pronounced as T is raised to 3.5 K. Again the ef- 
fects are most pronounced around tp ~ 7T. At higher 



temperatures, the supercurrent becomes too small to be 
detected with our method. We have obtained qualita- 
tively similar results on a shorter, but otherwise similar 
junction with A/e t h = 39. 

What is the origin of the microwave induced suppres- 
sion of the supercurrent around tp ~ 7r? The rf-field is 
expected to excite small oscillations of the magnetic flux 
enclosed by the loop [T3| ■ One possibility is that these os- 
cillations lead to the induction of an ac- voltage V ac across 
the junction and modify the distribution function /(e) of 
Andrccv states within the junction. One would then ex- 
pect sidebands in /(e) of width fiw r f [II]. Indeed, such a 
mechanism can lead to a distortion of the CPR qualita- 
tively similar to our observations. However, it disagrees 
with our data in that the I^\t) should again vanish 
exponentially for j > 1 and T > eth/ks- 

Another possibility are single-photon excitations of 
quasiparticles to ABS carrying supercurrent in the op- 
posite direction. Such excitations become possible when 
the minigap goes to zero near tp ~ n. To model this 
effect, we start, as in the equilibrium case, from Eq. [1] 
assuming that one can use the equilibrium supercurrent 
spectral density Imjs(e, tp) with a nonequilibrium distri- 
bution function /(e). 

A kinetic equation for /(e) in dirty superconductors 
under microwave irraduation was formulated first by 
Eliashberg [l9j . It can be represented in the form 

A{e) (/(e - u) - /(e)) + B(e) (/(e + u ) - /(e)) 
+ C(e)(l-f(u-e)-f(e))]=I[f(e)}, (3) 

where the coefficient P^ = D\5tpu/d e s\ 2 is proportional 
to the rf intensity. The phase oscillations Stp^ repre- 
sent the ac-response of the loop when excited by the mi- 
crowave field. The coefficients A(e), B(e), and C(e) can 
be expressed through the Green's functions of the sys- 
tem [19(. We have approximated the Green's functions of 
our proximity superconductor by the standard BCS form 
and replaced A with the minigap e g (tp). /[/(e)] is the 
scattering integral including electron-electron (ee) and 
electron-phonon (ep) interactions. The first two terms 
oc P u describe scattering of quasiparticles by photons and 
the third term pair-breaking, i.e., the creation of electron- 
hole pairs by single-photon absorption. 

Assuming that the ee-scattering is much stronger than 
the ep-scattering, we approximate /(e) by Fermi func- 
tions /*(e) with symmetrically shifted chemical poten- 
tials fj,* and —fi* [l0( for the microwave excited electron- 
and hole- like quasiparticles, respectively. Since the ex- 
citation of quasiparticles by photons creates no charge 
imbalance, the nonequilibrium distribution function /* 
is particle-hole symmetric. Then /* and fj,* can be de- 
termined from Eq. [3] within a relaxation time approxi- 
mation, rather than numerically solving the inhomogc- 
neous kinetic equation. We neglect the T-dependence of 



4 




a ^ =0) £/e Th b T(K) 



FIG. 4: (color) a) Spectral supercurrent js for tp/ir = 0.3, 0.6 
and 0.9 together with 1 — 2/* vs. e at 2.8 K. Inset: Zoom 
for e/e th < 3. Black lines: 1 - 2/* for P = 0, 24 dBm 
and 32 dBm (from top to bottom). The point 1 — 2/ = 
defines the chemical potential for electron-like quasi- 
particles. b) Normalized values of Iq (■,□), 7q 2 ' (#,0) and 
/q 3 ' (A, A) for the measured CPR without rf-irradiation (open 
symbols) and for 11.8 GHz and P = 28 dBm (filled symbols) 
as a function of temperature. The lines are the theoretical 
curves corresponding to f° and /*. 

the ep-scattering time r op in the limited T-interval un- 
der consideration. After integration of Eq. [3] over energy 
the terms containing A(e) and B(e) vanish, because they 
do not change the total number of excitations, and we 
obtain 

Pu> I™ C(e)(l-r(u-e)-r(e))de=— , (4) 

where An = J* °° (/*(e) - /°(e)) de and /°(e) is the equi- 
librium Fermi distribution function. 

In Fig. QJl we plot j s (e) together with 1 - 2/*(e). De- 
spite the resulting small values n* < 0.5 e t h, we have a 
noticable effect on the total supercurrent Isip) around 
ip ~ 7r. The increase of /i* cuts off the sharp positive 
peak of js at e s (ip) and enhances the relative weight of 
the negative part of js(e) at higher e in Eq. [TJ This 
can explain the observed reversal of Is{<p)- We express 
PuTph = n(ijj) exp (P[dBm]/10) in terms of a single pa- 
rameter rj, which also contains the w-dependent coupling 
between the sample loop and the antenna. P is the ex- 
ternally applied microwave power. The results of our 
calculations are displayed as solid lines in Fig. 3 with 
t]{uj) kept fixed to 0.001 for / = 7.0 GHz and 0.005 for 
11.8 GHz. All other parameters are identical to that of 
Fig. The agreement achieved is quite satisfactory. 

In Fig. [4}) we show the Fourier coefficients of the mea- 
sured CPR together with those extracted from the the- 
ory up to the 3 rd harmonics. For a better comparison, 
the data are normalized with respect to the equilibrium 
critical current Iq (T) . In thermodynamic equilibrium 
(open symbols), the sign of the Fourier coefficients alter- 
nates and they decay monotonically with temperature. 
Under rf-irradiation (full symbols), 1^ slightly decreases 



with microwave power, which may result from slight elec- 

(2) 

tron heating by the microwaves. In contrast, and 

(3) 

Jq are much larger, increase with T, and have opposite 
sign when compared to the case without rf irradiation. 
The experimental data agree rather well with the theo- 
retical values. Our observations resemble very much the 
peculiar temperature dependence of the amplitude of the 
subharmonic Shapiro steps observed in Refs. @, 01- 

In conclusion, we have experimentally verified long- 
standing theoretical predictions on the current-phase re- 
lation in diffusive Josephson junctions with highly trans- 
parent interfaces. We have presented strong evidence 
that microwave- irradiation of the junction induces higher 
harmonics in the current-phase relation. This effect can 
be understood in terms of a non-equilibrium occupa- 
tion of Andreev bound states carrying the supercurrent 
through the junction. Our results call for caution in the 
use of subharmonic Shapiro steps as a method for the 
determination of the current-phase relation. 
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